Symmetry breaking and self-interaction correction in the chromium atom and dimer.

Density functional approximations to the exchange-correlation energy can often identify strongly correlated systems and estimate their energetics through energy-minimizing symmetry-breaking. In particular, the binding energy curve of the strongly correlated chromium dimer is described qualitatively by the local spin density approximation (LSDA) and almost quantitatively by the Perdew-Burke-Ernzerhof generalized gradient approximation (PBE-GGA), where the symmetry breaking is antiferromagnetic for both. Here, we show that a full Perdew-Zunger self-interaction-correction (SIC) to LSDA seems to go too far by creating an unphysical symmetry-broken state, with effectively zero magnetic moment but non-zero spin density on each atom, which lies ∼4 eV below the antiferromagnetic solution. A similar symmetry-breaking, observed in the atom, better corresponds to the 3d↑↑4s↑3d↓↓4s↓ configuration than to the standard 3d↑↑↑↑↑4s↑. For this new solution, the total energy of the dimer at its observed bond length is higher than that of the separated atoms. These results can be regarded as qualitative evidence that the SIC needs to be scaled down in many-electron regions.

Orbital dependent complications for close vs well-separated electrons in diradicals.

We investigate two limits in open-shell diradical systems: O3, in which the interesting orbitals are in close proximity to one another, and (C21H13)2, where there is a significant spatial separation between the two orbitals. In accord with earlier calculations, we find that standard density-functional approximations do not predict the open-shell character for the former case but uniformly predict the open-shell character for the latter case. We trace the qualitatively incorrect behavior in O3 predicted by these standard density functional approximations to self-interaction error and use the Fermi-Löwdin-orbital-self-interaction-corrected formalism to determine accurate triplet, closed-shell singlet, and open-shell broken-spin-symmetry electronic configurations. Analysis of the resulting many-electron overlap matrices allows us to unambiguously show that the broken-spin-symmetry configurations do not participate in the representation of the Ms = 0 triplet states and allows us to reliably extract the singlet-triplet splitting in O3 by analyzing the energy as a function of Fermi-orbital-descriptor permutations. The results of these analyses predict the percentage of open-shell character in O3, which agrees well with conventional wavefunction-based methods. While these techniques are expected to be required in cases near the Coulson-Fischer point, we find that they will be less necessary in diradical systems with well-separated electrons, such as (C21H13)2. Results based on energies from self-interaction-corrected generalized gradient, local density, and Hartree-Fock approximations and experimental results are in generally good agreement for O3. These results help form the basis for deriving extended Heisenberg-like Hamiltonians that are needed for descriptions of molecular magnets when there are competing low-energy electronic configurations.

Use of FLOSIC for understanding anion-solvent interactions.

An Achille's heel of lower-rung density-functional approximations is that the highest-occupied-molecular-orbital energy levels of anions, known to be stable or metastable in nature, are often found to be positive in the worst case or above the lowest-unoccupied-molecular-orbital levels on neighboring complexes that are not expected to accept charge. A trianionic example, [Cr(C2O4)3]3-, is of interest for constraining models linking Cr isotope ratios in rock samples to oxygen levels in Earth's atmosphere over geological timescales. Here we describe how crowd sourcing can be used to carry out self-consistent Fermi-Löwdin-Orbital-Self-Interaction corrected calculations (FLOSIC) on this trianion in solution. The calculations give a physically correct description of the electronic structure of the trianion and water. In contrast, uncorrected local density approximation (LDA) calculations result in approximately half of the anion charge being transferred to the water bath due to the effects of self-interaction error. Use of group-theory and the intrinsic sparsity of the theory enables calculations roughly 125 times faster than our initial implementation in the large N limit reached here. By integrating charge density densities and Coulomb potentials over regions of space and analyzing core-level shifts of the Cr and O atoms as a function of position and functional, we unambiguously show that FLOSIC, relative to LDA, reverses incorrect solute-solvent charge transfer in the trianion-water complex. In comparison to other functionals investigated herein, including Hartree-Fock and the local density approximation, the FLOSIC Cr 1s eigenvalues provide the best agreement with experimental core ionization energies.

A Perspective on Sustainable Computational Chemistry Software Development and Integration.

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes.

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r2SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew-Burke-Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r2SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

Downward quantum learning from element 118: Automated generation of Fermi-Löwdin orbitals for all atoms.

A new algorithm based on a rigorous theorem and quantum data computationally mined from element 118 guarantees automated construction of initial Fermi-Löwdin-Orbital (FLO) starting points for all elements in the Periodic Table. It defines a means for constructing a small library of scalable FLOs for universal use in molecular and solid-state calculations. The method can be systematically improved for greater efficiency and for applications to excited states such as x-ray excitations and optically silent excitations. FLOs were introduced to recast the Perdew-Zunger self-interaction correction (PZSIC) into an explicit unitarily invariant form. The FLOs are generated from a set of N quasi-classical electron positions, referred to as Fermi-Orbital descriptors (FODs), and a set of N-orthonormal single-electron orbitals. FOD positions, when optimized, minimize the PZSIC total energy. However, creating sets of starting FODs that lead to a positive definite Fermi orbital overlap matrix has proven to be challenging for systems composed of open-shell atoms and ions. The proof herein guarantees the existence of a FLOSIC solution and further guarantees that if a solution for N electrons is found, it can be used to generate a minimum of N - 1 and a maximum of 2N - 2 initial starting points for systems composed of a smaller number of electrons. Applications to heavy and super-heavy atoms are presented. All starting solutions reported here were obtained from a solution for element 118, Oganesson.

Density Matrix Implementation of the Fermi-Löwdin Orbital Self-Interaction Correction Method.

The Fermi-Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi-Löwdin orbitals which are used to remove the self-interaction error in the Perdew-Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi-Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn-Sham with one additional term in the Kohn-Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn-Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. We analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville-von Neumann equation to calculate excitation energies.

Complex Fermi-Löwdin orbital self-interaction correction.

This paper introduces the use of complex Fermi orbital descriptors (FODs) in the Fermi-Löwdin self-interaction-corrected density functional theory (FLOSIC). With complex FODs, the Fermi-Löwdin orbitals (FLOs) that are used to evaluate the SIC correction to the total energy become complex. Complex FLO-SIC (cFLOSIC) calculations based on the local spin density approximation produce total energies that are generally lower than the corresponding energies found with FLOSIC restricted to real orbitals (rFLOSIC). The cFLOSIC results are qualitatively similar to earlier Perdew-Zunger SIC (PZ-SIC) calculations using complex orbitals [J. Chem. Phys. 80, 1972 (1984); Phys. Rev. A 84, 050501(R) (2011); and J. Chem. Phys. 137, 124102 (2012)]. The energy lowering stems from the exchange-correlation part of the self-interaction correction. The Hartree part of the correction is more negative in rFLOSIC. The energy difference between real and complex solutions is greater for more strongly hybridized FLOs in atoms and for FLOs corresponding to double and triple bonds in molecules. The case of N2 is examined in detail to show the differences between the real and complex FLOs. We show that the complex triple-bond orbitals are simple, and physically appealing combinations of π and σg orbitals that have not been discussed before. Consideration of complex FODs, and resulting unitary transformations, underscores the fact that FLO centroids are not necessarily good guesses for FOD positions in a FLOSIC calculation.

Self-interaction correction in water-ion clusters.

We study the importance of self-interaction errors in density functional approximations for various water-ion clusters. We have employed the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in conjunction with the local spin-density approximation, Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and strongly constrained and appropriately normed (SCAN) meta-GGA to describe binding energies of hydrogen-bonded water-ion clusters, i.e., water-hydronium, water-hydroxide, water-halide, and non-hydrogen-bonded water-alkali clusters. In the hydrogen-bonded water-ion clusters, the building blocks are linked by hydrogen atoms, although the links are much stronger and longer-ranged than the normal hydrogen bonds between water molecules because the monopole on the ion interacts with both permanent and induced dipoles on the water molecules. We find that self-interaction errors overbind the hydrogen-bonded water-ion clusters and that FLOSIC reduces the error and brings the binding energies into closer agreement with higher-level calculations. The non-hydrogen-bonded water-alkali clusters are not significantly affected by self-interaction errors. Self-interaction corrected PBE predicts the lowest mean unsigned error in binding energies (≤50 meV/H2O) for hydrogen-bonded water-ion clusters. Self-interaction errors are also largely dependent on the cluster size, and FLOSIC does not accurately capture the subtle variation in all clusters, indicating the need for further refinement.

Electromagnetic control of spin ordered Mn3 qubits: a density functional study.

[Mn3O(O2CMe)(dpd3/2)]2 is composed of two monomers each of which contain three Mn atoms at the vertices of an equilateral triangle. A full analysis of the electronic and magnetic structure of the dimer shows that each Mn atom carries a local spin of S = 2 while other spin states are energetically much higher. This result suggests application for conventional as well as quantum tasks. A detailed analysis of the electronic and magnetic structure of the monomer, on the other hand, suggests that there are three spin states of S = 1, S = 3/2 and S = 2 per monomer which are energetically competitive. We found that while monomer-monomer interactions are very weak, the coupling of monomers via covalent linkers affects both the magnetization and electronic energy levels of monomers. In particular, the isolated monomers prefer a ground state with local spin of S = 1 on Mn atoms and an antiferromagnetically ordered structure while the dimers possess a ground state with local spin of S = 2 on Mn atoms and a ferromagnetically ordered structure. The investigation of the polarizability of both monomer and dimer is examined for antiferromagnetically ordered structures which induces a high dipole moment of 0.08 (a.u.) and 0.16 (a.u.) for monomer and dimer, respectively. The energy of the antiferromagnetic structure is also high compared to other spin-electric molecules.

A multiferroic molecular magnetic qubit.

The chiral Fe3O(NC5H5)3(O2CC6H5)6 molecular cation, with C3 symmetry, is composed of three six-fold coordinated spin-carrying Fe3+ cations that form a perfect equilateral triangle. Experimental reports demonstrating the spin-electric effect in this system also identify the presence of a magnetic uniaxis and suggest that this molecule may be a good candidate for an externally controllable molecular qubit. Here, we demonstrate, using standard density-functional methods, that the spin-electric behavior of this molecule could be even more interesting as there are energetically competitive reference states associated with both high and low local spins (S = 5/2 vs S = 1/2) on the Fe3+ ions. Each of these structures allow for spin-electric ground states. We find that qualitative differences in the broadening of the Fe(2s) and O(1s) core levels, shifts in the core-level energies, and the magnetic signatures of the single-spin anisotropy Hamiltonian may be used to confirm whether a transition between a high-spin manifold and a low spin manifold occurs.

Magnetic Signatures of Hydroxyl- and Water-Terminated Neutral and Tetra-Anionic Mn12 -Acetate.

We investigate the energetics and magnetic signatures of the parent molecular magnet Mn12 -Acetate [Mn12 O12 (COOR)16 (H2 O)4 ] and a chemically decomposed version of this structure, in which the four water molecules are converted to hydroxyl groups and hydrogen molecules. We determine electron addition and water decomposition energetics for this water-containing molecule using density-functional methods and include the recent Fermi-Löwdin-Orbital self-interaction correction. We find that it only costs 0.32 eV to add four electrons to the parent molecule. Furthermore, due to the strong Coulomb attractions between hydroxyl anions and the Mn cations, the energy cost for breaking the four coordinating water molecules into four coordinating hydroxyls and two hydrogen molecules is decreased in the tetra-anionic structure relative to the neutral structure. We calculate magnetic anisotropy barriers for the neutral molecule and the dehydrogenated tetra-anion and show that large changes in the magnetic anisotropy arise the strong attraction between the hydroxyl anions and four of the crown Mn cations. We suggest that the large changes in magnetic signals associated with the [Mn12 O12 (COOR)16 (HOH)4 ] to [Mn12 O12 (COOR)16 (OH- )4 + 2H2 ] decomposition could provide a nondestructive spectroscopic method for learning about water decomposition mechanisms in a class of realizable model catalytic systems that have been synthesized recently. © 2019 Wiley Periodicals, Inc.

Theoretical studies of the vibrational properties of octahedrane (C12H12): A polyhedral caged hydrocarbon molecule.

The vibrational properties of octahedrane (C12H12) are calculated using density-functional theory employing two different computational methods: an all-electron Gaussian orbital approach and a Naval Research Laboratory-tight-binding scheme (NRL-TB) coupled with molecular dynamics (NRL-TBMD). Both approaches yield vibrational densities of states for octahedrane that are in good general agreement with each other. NRL Molecular Orbital Library can also provide accurate infrared and Raman spectra which can be analyzed and compared with experimental results, while NRL-TBMD can be conveniently scaled up for larger finite-temperature simulations. This latter approach is used in our paper to produce a theoretical prediction for a stable room temperature structure of octahedrane.

Symmetry Breaking within Fermi-Löwdin Orbital Self-Interaction Corrected Density Functional Theory.

Fermi-Löwdin orbital self-interaction corrected density functional theory (FLO-SIC DFT) is applied to C3H5, NO3-, O3, and CH3. In general our results indicate that FLO-SIC does favor symmetric setups for molecules with nontrivial chemical bonding. Further we discuss two types of possible symmetry breaking. In the case of O3 a ground state density that breaks symmetry can be found for the symmetric molecular geometry that may be caused by insufficient treatment of correlation energy. The CH3 radical presents a second type of symmetry breaking were the lowest energy geometry becomes distorted. The latter highlights the importance of further development of approximate DFT functionals as well as further extensions of the FLO-SIC method to overcome such nonphysical artifacts.

Self-interaction corrections applied to Mg-porphyrin, C60, and pentacene molecules.

We have applied a recently developed method to incorporate the self-interaction correction through Fermi orbitals to Mg-porphyrin, C60, and pentacene molecules. The Fermi-Löwdin orbitals are localized and unitarily invariant to the Kohn-Sham orbitals from which they are constructed. The self-interaction-corrected energy is obtained variationally leading to an optimum set of Fermi-Löwdin orbitals (orthonormalized Fermi orbitals) that gives the minimum energy. A Fermi orbital, by definition, is dependent on a certain point which is referred to as the descriptor position. The degree to which the initial choice of descriptor positions influences the variational approach to the minimum and the complexity of the energy landscape as a function of Fermi-orbital descriptors is examined in detail for Mg-porphyrin. The applications presented here also demonstrate that the method can be applied to larger molecular systems containing a few hundred electrons. The atomization energy of the C60 molecule within the Fermi-Löwdin-orbital self-interaction-correction approach is significantly improved compared to local density approximation in the Perdew-Wang 92 functional and generalized gradient approximation of Perdew-Burke-Ernzerhof functionals. The eigenvalues of the highest occupied molecular orbitals show qualitative improvement.

Fermi orbital self-interaction corrected electronic structure of molecules beyond local density approximation.

The correction of the self-interaction error that is inherent to all standard density functional theory calculations is an object of increasing interest. In this article, we apply the very recently developed Fermi-orbital based approach for the self-interaction correction [M. R. Pederson et al., J. Chem. Phys. 140, 121103 (2014) and M. R. Pederson, J. Chem. Phys. 142, 064112 (2015)] to a set of different molecular systems. Our study covers systems ranging from simple diatomic to large organic molecules. We focus our analysis on the direct estimation of the ionization potential from orbital eigenvalues. Further, we show that the Fermi orbital positions in structurally similar molecules appear to be transferable.

Communication: practical and rigorous reduction of the many-electron quantum mechanical Coulomb problem to O(N(2/3)) storage.

It is tacitly accepted that, for practical basis sets consisting of N functions, solution of the two-electron Coulomb problem in quantum mechanics requires storage of O(N(4)) integrals in the small N limit. For localized functions, in the large N limit, or for planewaves, due to closure, the storage can be reduced to O(N(2)) integrals. Here, it is shown that the storage can be further reduced to O(N(2/3)) for separable basis functions. A practical algorithm, that uses standard one-dimensional Gaussian-quadrature sums, is demonstrated. The resulting algorithm allows for the simultaneous storage, or fast reconstruction, of any two-electron Coulomb integral required for a many-electron calculation on processors with limited memory and disk space. For example, for calculations involving a basis of 9171 planewaves, the memory required to effectively store all Coulomb integrals decreases from 2.8 Gbytes to less than 2.4 Mbytes.

Fermi orbital derivatives in self-interaction corrected density functional theory: Applications to closed shell atoms.

A recent modification of the Perdew-Zunger self-interaction-correction to the density-functional formalism has provided a framework for explicitly restoring unitary invariance to the expression for the total energy. The formalism depends upon construction of Löwdin orthonormalized Fermi-orbitals which parametrically depend on variational quasi-classical electronic positions. Derivatives of these quasi-classical electronic positions, required for efficient minimization of the self-interaction corrected energy, are derived and tested, here, on atoms. Total energies and ionization energies in closed-shell singlet atoms, where correlation is less important, using the Perdew-Wang 1992 Local Density Approximation (PW92) functional, are in good agreement with experiment and non-relativistic quantum-Monte-Carlo results albeit slightly too low.

Communication: self-interaction correction with unitary invariance in density functional theory.

Standard spin-density functionals for the exchange-correlation energy of a many-electron ground state make serious self-interaction errors which can be corrected by the Perdew-Zunger self-interaction correction (SIC). We propose a size-extensive construction of SIC orbitals which, unlike earlier constructions, makes SIC computationally efficient, and a true spin-density functional. The SIC orbitals are constructed from a unitary transformation that is explicitly dependent on the non-interacting one-particle density matrix. When this SIC is applied to the local spin-density approximation, improvements are found for the atomization energies of molecules.

Photoelectron spectroscopic and computational studies of the Pt@Pb10⁻¹ and Pt@Pb12¹â»/²â» anions.

A combination of anion photoelectron spectroscopy and density functional theory calculations has elucidated the geometric and electronic structure of gas-phase endohedral Pt/Pb cage cluster anions. The anions, Pt@Pb10⁻¹ and Pt@Pb12¹â» were prepared from "preassembled" clusters generated from crystalline samples of [K(2,2,2-crypt)]2[Pt@Pb12] that were brought into the gas phase using a unique infrared desorption/photoemission anion source. The use of crystalline [K(2,2,2-crypt)]2[Pt@Pb12] also provided access to K[Pt@Pb(n)](-) anions in the gas phase (i.e., the K⁺ salts of the Pt@Pb(n)²â» anions). Anion photoelectron spectra of Pt@Pb10⁻¹, Pt@Pb12¹â», and K[Pt@Pb12]¹â» are presented. Extensive density functional theory calculations on Pt@Pb10³â»/²â»/¹â»/° and Pt@Pb12²â»/¹â» provided candidate structures and anion photoelectron spectra for Pt@Pb10⁻¹ and Pt@Pb12¹â». Together, the calculated and measured photoelectron spectra show that Pt@Pb10⁻¹ and Pt@Pb12²â»/¹â» endohedral complexes maintain their respective D(4d) and slightly distorted I(h) symmetries in the gas phase even for the charge states with open shell character. Aside from the fullerenes, the Pt@Pb12²â» endohedral complex is the only bare cluster that has been structurally characterized in the solid state, solution, and the gas phase.